Differential Dynamic Microscopy of Bacterial Motility 
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We demonstrate 'differential dynamic microscopy' (DDM) for the fast, high throughput char- 
acterization of the dynamics of active particles. Specifically, we characterize the swimming speed 
distribution and the fraction of motile cells in suspensions of Escherichia coli bacteria. By averaging 
over ~ 10 4 cells, our results are highly accurate compared to conventional tracking. The diffusivity 
of non-motile cells is enhanced by an amount proportional to the concentration of motile cells. 



Diverse processes in multicellular organisms such as 
chemotaxis involve motility which is also ubiquitous 
in unicellular organisms such as bacteria, enabling, e.g. 
the pathogen Helicobacter pylori to invade the stomach 
epithelium 2]. Globally, bacterial motility may be cou- 
pled to aquatic nutrient recycling Q. The bacterium 
Escherichia coli is a paradigm for understanding cell 
motility [J. A cell executes a random walk by alternat- 
ing between swimming (or 'running') at average speed 
v > 10/im/s for ~ 1 s and tumbling for ~ 0.1 s. 

Early bacterial motility work relied on tracking one 
to a few cells 0, S]. Today, ~ 10 2 - 10 3 cells can be 
tracked simultaneously Tracking yields a host of 

parameters, including v (e.g. in polymer solutions @) 
and the fraction of motile organisms, a (e.g. in oceanic 
bacteria Q). But tracking is laborious, and the need for 
averaging over many data sets to achieve high accuracy 
restricts the scope for time-dependent measurements. 

We demonstrate a fast, high throughput method for 
characterizing E. coli motility. It should be applicable to 
other bacteria and micro-organisms, and to a new gener- 
ation of synthetic, self-propelled 'active particles' (lfj{ . 

Dynamic light scattering (DLS), long used for measur- 
ing diffusivity in colloids, is in principle suitable for the 
fast characterization of motile bacteria [ll|. DLS yields 
the normalized intermediate scattering function (ISF), 
f(q, r) (where q is the scattering vector and r is time) 
[12j, which probes density relaxation processes at length 
scale 27r/g. But the lowest scattering angle in conven- 
tional DLS, ~ 20° (or q ~ 4.5^m _1 ), probes dynamics 
at 27r/g < 1.4/im, where cell body precession [13| and 
other motions in E. coli contribute strongly to the decay 
of the ISF. Thus, contrary to initial claims [ll|, E. coli 
swimming, which occurs on the scale of v/r Iun ~ 10pm, 
cannot be characterized unambiguously using DLS unless 
we can access q < 27r/10pm ~ O.Bpm^ 1 (or < 3°) [l3j [. 

Instead of implementing such ultra- low-angle DLS, we 
use the powerful technique of Differential Dynamic Mi- 
croscopy (DDM) to measure f(q,r) for bacterial swim- 
ming. A form of DDM was first used to study density 
fluctuations in binary mixtures |14j . It has recently been 
used to measure colloidal diffusivity [15j], requiring only 
non-specialized equipment (microscope, camera and com- 



puter). The DDM of colloids, however, does not utilize its 
unique capability to reach very low q (< lpm -1 ), which 
turns out to be essential for probing bacterial swimming. 

The theory of DDM is detailed in We give an 

alternative derivation, which also explains experimental 
procedures. The raw data are time-lapsed images of (say) 
bacteria, described by the intensity I(r, t) in the image 
plane (f). From these we calculate difference images at 
various delay times, r, D{r,r) = I(f,t + r) — I{r,t) = 
AI(r, t) - AI(r, 0), where AI(r, t) = I(r, t) - {I) denotes 
intensity fluctuations. Fourier transforming D(r, r) gives 



F D {q,T) = J D{f,T)e l ^dr. 



(1) 



For stationary, isotropic processes, we average over the 
start time t in the difference images and azimuthally in 
q space to calculate the basic output of DDM, what we 
may call the 'differential intensity correlation function' 
(DICF), (\F D (q, r)| 2 ) (where q = \q\). 

We now show that the DICF is related simply to the 
ISF if we assume that intensity fluctuations in the image 
are proportional to the fluctuations in the number density 
of bacteria around the average density (p) : 



AI(f,t) = K,Ap(r,t) 



(2) 



Here the constant k depends on the contrast mechanism 
and Ap(f, t) = p(r, t) — (p). Eqs. (P) and © now give 



F D (q,r) = K[Ap(q,T) - Ap(q,0)} 



where Ap{q, r) 



Ap(f, t)e l ^&r 



Thus, the DICF can be expressed as 

(Ap(q,0)Ap(q,r)) 



(\F D (q,r)\ 2 ) = A(q) 



1 - 



[Ap(?)] : 



where A(q) = 2K 2 ([Ap(q)]' 



(3) 
(4) 

(5) 
(6) 



The prefactor A{q) depends on the imaging system, k, 
and on the sample's structure, ([Ap(q)} 2 }. Recogniz- 
ing that the r-dependent term on the right-hand side 
of Eq. ([S]) is the ISF, we arrive at this key result: 



(\F D (q,r)\ 2 ) = A(q) [1 - f(q, r)] + B(q) 
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FIG. 1: Reconstructed ISFs of non-motile bacteria plotted 
against q 2 r. Solid (black) curves for over 200 values of q in 
the range 0.5/^m -1 < q < 2.2/im -1 collapse, but curves from 
lower q (grey) do not collapse. Inset: fitted diffusivity D(q) 
(black and grey with the same meaning). 



where we have included a term B{q) to account for cam- 
era noise. Thus, the power spectrum of intensity fluc- 
tuations of the images, (\Fjj{q, r)| 2 ), yields the ISF. In 
practice, we reconstruct f{q,T) by using a parametrized 
model of the ISF to fit the measured DICF with Eq. Q. 

For independent particles, f{q,r) = ( e ~ I 9' AT '(' r )), 
where Ar(r) is the single-particle displacement 12]. This 
reduces to f{q, t) = e~ Dq T for identical diffusing spheres 
with diffusivity D [121] . For a swimmer with velocity v, 
Ar(r) = vt. For an isotropic population of such swim- 
mers in 3D, f(q,r) = sm{qvr)/qvT = sinc(gi)r) 
since a swimmer inevitably also undergoes Brownian mo- 
tion, this needs to be multiplied by an exponential pref- 
actor e~ Dq ~ T . If only a fraction a of swimmers are motile 
with speed distribution P{v), then the full ISF reads (l7j : 



/(?.r) 



-Dq T 



:i 



P(v)sinc(qvT)dv 



(8) 



In order to use this model to interpret our DDM data 
from E. coli, we need to specify a form for P(v). Limited 
previous data [TJ [I?} suggest a peaked function with 
P(v — » 0) — >■ 0. We use a Schulz distribution: 
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v 



1) 



(9) 



where Z is related to the variance a 2 of the distribution 
by a = v{Z + l) -1 / 2 . The integral in Eq. ([8]) evaluates 
toll 



P{v)smc{qvr)d 



sin ( Z tan 



where = (qvr)/(Z + 1). 



-(10) 
(11) 



We studied E. coli AB1157 grown at 30°C in L- 
broth, re-inoculated into T-broth and harvested in mid- 
exponential phase, washed three times by filtration 



FIG. 2: Symbols: reconstructed ISFs for motile E. coli at 
8 values of q spanning the whole available range (see key). 
Lines: calculated ISFs, Eqs. (I5)l- (llip . using fitted parameters. 



(0.45/im filter) in motility buffer and re-suspended in the 
same buffer to an optical density of 0.3 (at 600nm), giv- 
ing a final cell volume fraction of <f> w 0.06%. (See Supple- 
mentary Material for details.) Care was taken through- 
out to minimize damage to fiagella. A ~ 400/^m deep 
flat glass cell was filled with ~ 150/d of cell suspension, 
sealed, and observed at 22±1°C. Swimming behavior was 
constant over a 15 minute period. We also used a non- 
motile mutant with 'paralyzed' fiagella {mot A). 

We collected movies of cells using a 10 x phase- 
contrast objective in a Nikon Eclipse Ti inverted micro- 
scope. Images were obtained « 100/im from the bot- 
tom of a 400/im-thick sample cell. A high-speed cam- 
era (Mikrotron MC 1362) was connected to a PC with a 
frame grabber card with 1GB onboard memory. Movies 
were acquired typically at 100Hz. The frame size L 2 
was 500 x 500 and 1024 x 1024 pixels for motA mutants 
and wild-type cells respectively, imaging ~ 10 4 cells in a 
0.7mm 2 or 1.4mm 2 field of view over 38s or 8s. The pixel 
size (or spatial sampling frequency) is k — 0.712/zm -1 , 
so that <7 m i n = 2-nk/L ss O.Ol^m -1 or 0.004 ^m _1 . 

To calculate the DICFs from the raw images, we used 
a Lab View (National Instruments) code optimized for 
an 8-core PC (dual Intel Xeon quad-core processors, 
2GHz/core, 4GB RAM). Analyzing ~ 40s of movies takes 
~ 10 minutes. We then fitted each DICF to Eq. ([7]) using 
Eqs. (I5|- (|lip . At each q, non-linear least-squares fitting 
using the Levenberg-Marquardt algorithm (l9| in IGOR 
Pro (WaveMetrics) returns six parameters: v, a, D, a, 
A and B. Fitting the whole q range takes ~ 30 s. From 
the fitted A{q) and B{q), we obtain the reconstructed 
ISF using the measured DICF and Eq. ([7]). We also ob- 
tain the calculated ISF by using the fitted {v, a, D, a} in 
Eqs. ©-(UTJ. 

We first studied non-motile {motA) cells. Measured 
DICFs are well fitted using Eq. © with f(q, r) = er D ^ T 
(i.e. Eq. © with a = 0) (Supp. Fig. 1). The fitted 
diffusivity, D{q), was q independent within experimen- 
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FIG. 3: Parameters extracted from fitting the DICF data for 
motile cells. From top to bottom: v and er of the Schulz 
distribution, motile fraction a and diffusivity D. 



tal uncertainties in the range 0.5/im _1 < q < 2.2/zm _1 , 
Fig. [T] inset, and averaged to _D = 0.30 ± 0.01/xm 2 /s. 
Conventional DLS (data not shown) gave an exponential 
f(q, t) and D = 0.32 ± 0.02^m 2 /s, agreeing with DDM. 
The reconstructed ISFs collapse onto each other in the 



range 0.5/mi~ 



1 < q < 2.2/inr 



when plotted against q r 



(black curves, Fig. Q}, i.e. the non-motile cells are purely 
diffusive. At q < 0.5/im _1 , f(q,r) has not decayed to 
zero at the longest time probed in our experiments, so 
that fitting becomes less reliable because of the difficulty 
in estimating A(q) (cf. Eq. (O). The reconstructed ISFs 
therefore do not collapse under g 2 r scaling and D(q) is 
noisy (grey curves, Fig.[TJ crosses, Fig. Q] inset). 

We next studied motile cells. The measured DICFs 
(Supp. Fig. 2) were again fitted to Eq. ([7]), now using 
the full f(q,r) in Eq. © and a Schulz P(v), Eqs. ©- 
(fTTjl . A selection of the reconstructed ISFs is shown in 
Fig. [U (points), where we also superimpose the calculated 
ISFs (curves). The ISFs display a characteristic shape, 
especially at low q: a fast decay dominated by swim- 
ming followed by a slower decay dominated by diffusion. 
(Compare also the different time axes in Figs. [T]and[2J) 

All fit parameters characterizing swimming are shown 
in Fig. [3] [25|. The noise increases at low q, primarily 
because the long-time, diffusive part of /((/,t) has not 
reached zero in our time window at these q, Fig. [3J ren- 
dering it harder to determine the diffusivity accurately: 
the \ow-q noise is particularly evident in the fitted D(q), 
Fig. [3J But to within experimental uncertainties all pa- 
rameters in Fig. P are essentially ^-independent at least 



for 



> 



1/im 



suggesting that our model is able 



to capture essential aspects of the dynamics of a mixed 
population of non-motile and motile E. coli. Averaging 
yields v = 13.7 ± 0.1/xms _1 and a — 7.0 ± 0.1/zms -1 , 
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FIG. 4: ■: D as a function of a (lower axis) and the volume 
fraction of swimmers, <f) a (upper axis). The black line is a 
linear fit. •: D(a) from simulations. 



with error bars reflecting estimated residual q dependen- 
cies. Changing A(q) and B{q) by using a 20 x objective 
(which is sub-optimal for our experiment) produced the 
same fitted motility parameters in the relevant q range. 

The presentation so far rests on a number of assump- 
tions. Our derivation of Eq. ([7]) assumes that the decor- 
relation of f(q,r) caused by the change in intensity of a 
swimmer's image due to its motion along the optic (z) 
axis can be neglected. While wild-type E. coli AB1157 
tumbles between 'runs' and the swim path between tum- 
bles is slightly curved, Eq. (J8j) neglects these effects. We 
tested the validity of these assumptions by performing 
DDM analysis on the output from computer simulations. 

We carried out Brownian dynamics simulations of non- 
interacting point particles at a number density and in 
a geometry directly comparable to our experiments. A 
fraction a of the particles had a drift speed drawn from 
a Schulz distribution. From these simulations, we con- 
structed a sequence of 2D pixellatcd 'images' with the 
same field of view as in experiments. All particles in a 
slice of thickness d centered at z = contribute to the 
image. A particle at {x, y, z) is 'smeared' into an 'image' 
covering the pixel containing (x, y) and its 8 neighboring 
pixels. The contrast of the image, c, depends on z. We 
experimentally determined d and c(z)by imaging a single 
bacterium as the focal plane traversed the cell. The mea- 
sured c(z) could be fitted by a symmetric quadratic that 
dropped to background noise outside a ~ 40/mi slice. 

As input, we used v = 13.7/xms^ 1 , a = 7.0/ims _1 , 
a = 0.577 and D = 0.543/mi 2 /s (cf. Fig. [3]). Fit- 
ting DICFs calculated from simulated 'images' (Supp. 
Fig. 3a) gave g-independent outputs (Supp. Fig. 3b): 
v_ = 13.8 ± 0.1, o = 7.2 zb 0.2, a = 0.58 ± 0.01 and 
D = 0.55 ± 0.02 (where the uncertainties are standard 
deviations), agreeing with inputs. Thus, at d = 40/im 
depth of field, the intensity decorrelation due to z mo- 
tion has negligible effect, presumably because it is much 
slower than the decorrelation due to swimming and diffu- 
sion. However, if we scale c{z) to smaller depths of field, 
the fitting eventually fails at d « 10/mi (data not shown): 
at this small focus depth, a small z movement produces 
a large intensity variation, invalidating our analysis. 
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DDM essentially determines the (inverse) time it takes 
a cell to traverse ~ 2n/q, i.e. it measures 'linear speeds'. 
Tumbling or curvature will therefore lower the measured 
speed, with the effect more noticeable at lower q. Our 
experimental v(q), Fig. [3] does indeed show a slight de- 
crease towards low q. As expected, however, the v{q) 
recovered from analysing simulated straight swimmers 
(Supp. Fig. 3) show no such dependence. More de- 
tailed analysis of the measured v(q) may therefore yield 
further information about tumbling and curvature. 

We next mixed suspensions of bacteria with known 
a with various proportions of non-motile cells, creating 
samples with < a < 0.8. DDM shows that D in- 
creases with a, Fig. |4] (squares). Since the fitting of D 
from Eq. [8] is largely determined by the diffusion of the 
non-swimmers, Fig. U shows that swimmers enhance the 
diffusion of non-swimmers. Note that this observation is 
not a fitting artifact: such dependence is not observed in 
simulations, which returned an a-independcnt D, Fig. 2] 
Since the simulated particles are non-interacting, our ex- 
perimental observation must be due to direct or hydrody- 
namic interaction between swimmers and non-swimmers. 

The enhanced diffusion of (passive) particles in sus- 
pensions of motile E. coli has been observed before using 
direct tracking at both low concentration (tp = 0.003%) 
in 3D 20] and high concentration (0 « 10%) in 2D 
Scaling arguments suggest that in the limit of indepen- 
dent swimmers, the enhancement should scale linearly as 
the concentration of swimmers (f> a (22J. In our experi- 
ments, the volume fraction of non-swimmers varies, but 
remains < 0.1%, i.e., they can be considered as indepen- 
dent 'tracer' particles. Thus, Fig. [4] can be reinterpreted 
as a plot of the effective diffusion coefficient of tracer 
particles as the concentration of swimmers increases from 
<j> a =Qto<j) a = 0.06% x 0.8 = 0.048%. The enhancement 
indeed scales linearly with the swimmer concentration. 



Our D{a) results may also be compared to enhanced 
tracer diffusion by Chlamydomonas reinhardtii, a nearly- 
spherical single-cell algae [23T ], However, C. reinhardtii 
(radius ~ 5/im) are larger than E. coli, and swim much 
faster (y ~ 100 /urns -1 ). More fundamentally, C. rein- 
hardtii (a 'puller') and E. coli (a 'pusher') generate qual- 
itatively different flow fields, which may have conse- 
quences for tracer diffusion [22J. Nevertheless, it is in- 
triguing that 2% of C. reinhardtii quadruples the dif- 
fusivity of 2/mi tracers, while 0.048% of motile E. coli 
already doubles the diffusivity of non-motile cells. 

To summarize, we have shown that DDM is a fast, 
high-throughput method for characterizing the bulk 
motility of wild-type E. coli. The method could, in prin- 
ciple, be extended to characterize the run-tumble-run 
random walk of individual cells [24} (by going to even 
lower q), or to the study of motility near surfaces (which 
requires the use of a different f(q, r) in Eq. ©). The 
method may also be applicable to the study of other 
motile organisms, including spermatozoa, as well as for 
characterizing the motions of synthetic motile colloids 
[lcj . But the q range, camera speed and data acquisition 
time need to be optimized for each particular class of 
motility to be characterized. Our finding that even low 
concentrations of motile cells enhance the diffusivity of 
non-motile cells may have implications for understanding 
the coupling between bacterial motility and the recycling 
of organic debris in natural aqueous habitats Q . 

Finally, we should emphasize that DDM yields f(q,r) 
of suspensions of active swimmers irrespective of <fi, pro- 
vided that Eq. @ remains valid. It is therefore a general 
method for studying the dynamics of these suspensions, 
including interaction effects at higher <f>, although new 
models will clearly be needed for interpreting the data. 
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